Thermal Fluctuations of Vortex Matter in Trapped Bose— Einstein Condensates 
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We perform Monte Carlo studies of vortices in three dimensions in a cylindrical confinement, with 
uniform and nonuniform density. The former is relevant to rotating ^He, the latter is relevant to a 
rotating trapped Bose-Einstein condensate. In the former case we find dominant angular thermal 
vortex fluctuations close to the cylinder wall. For the latter case, a novel effect is that at low 
temperatures the vortex solid close to the center of the trap crosses directly over to a tension-less 
vortex tangle near the edge of the trap. At higher temperatures an intermediate tensionful vortex 
liquid located between the vortex solid and the vortex tangle, may exist. 



Topological vortex excitations are hallmarks of quan- 
tum ordered states such as superconductors and super- 
fiuids. Vortices are important at vastly different length 
scales ranging from the dynamics of neutron stars to 
transport properties in superconductors and rotational 
response of ultracold atomic gases in optical traps |l|, [3 • 
Particularly rich physics is associated with various or- 
der ings of quantum vortices, and the transitions between 
them. In certain systems, even numerous aggregate 
states of vortex matter are possible 3]. The simplest, yet 
extremely important, temperature-induced phase transi- 
tion in vortex matter in a one-component superconductor 
or superfluid is the one from an ordered vortex line lat- 
tice (VLL) to a disordered vortex liquid (VL) state. This 
has been a subject of intense research in the context of 
high- Tc superconductivity. 

Remarkable progress has recently been made in the 
physics of vortices in Bose-Einstein condensates (BEC) 
of ultra-cold atoms jjj, I3 • These systems are extremely 
clean and have physical parameters controllable in a wide 
range. This has led to many suggestions for testing a 
number of general physical concepts in BECs. Since in- 
dividual vortices, as well as their ordering patterns, can 
be resolved in ultra-cold gases HI 5 a natural question 
arises if in these systems a thermally induced phase tran- 
sition from VLL to VL can be observed, and if it can shed 
more light on fundamental properties of the VL states in 
quantum fluids in general. Detailed theoretical studies 
of a model uniform system, with density corresponding 
to the values attainable in the center of a trap, indicate 
that in present experiments rotation is too slow and/or 
particle numbers are too high to obtain a VLL melting 
if the approximation of a uniform system holds |4^ . We 
also mention studies of quantum VLL melting in a one 
dimensional optical lattice p]. In a harmonic trap, the 
condensate density will gradually be depleted from the 
center of the trap towards its edge. This suggests that 
thermal fluctuation effects will be enhanced close to the 
edges. We thus focus on investigating possible crossover 
states of vortex matter in a spatially inhomogeneous sys- 
tem such as a trapped BEC. 

Melting of vortex lines in a bosonic condensate can 



be modeled by the uniformly frustrated lattice 3D XY 
model (see e.g. 0,13,13). We study a system which takes 
into account the presence of a harmonic trap by using 
the Hamiltonian 



H 



(ij) 



cos{Oj — Oi — Aij), 



(1) 



where Oi is the phase of the condensate wave function at 
position i. The factor Pij = P{rij) = 1 — (r^j/i?)^ if Vij < 
R^ where R is the cloud size and rij is the radial distance 



from the cloud center. For 



> R we set Pij = 



0. 



We account for circulation by introducing the potential 
Aij = J- dl- A. Here, V x A = (0, 0, 27r/), and / is the 
number of rotation-induced vortices per plaquette in the 
x7/-plane. 

We perform Monte Carlo (MC) simulations on Eq. (^ 
with the Metropolis algorithm. We have considered cu- 
bic numerical grids of size L^, with L = 36, 72, and along 
the z-direction we impose periodic boundary conditions 
to model an elongated system. The filling fraction is 
/ = 1/36, and temperatures T G (0.30,3.0) 0. The 
temperature at which the Bose condensation takes place 
at zero rotation in a bulk system, is given by T = 2.2 
in our units. We note that a rapidly rotating system 
should be better described by the Lowest Landau Level 
(LLL) approximation 10], characterized for instance by 
the fact that the intervortex separation is comparable to 
the vortex-core size. However, several experiments on 
elongated systems are clearly outside the LLL regime. 
For instance, in Ref. 2], the inter- vortex separation is 
given as 5.0/i7Ti, while the healing length (vortex-core ra- 
dius ^) is 0.2/i7Ti. A parameter estimating the validity of 
the LLL approximation is the ratio A of the interaction 
energy scale to the level spacing of the transverse har- 
monic confinement 10], A = 47rh'^asn/{Mhu±)^ where a^ 
is the s-wave scattering length, M is the particle mass, 
and u;± is the trap frequency. The LLL approximation 
requires A <C 1, while the parameters of Ref. 2] corre- 
spond to at least A G 1 — 100. Under such conditions, we 
expect Eq. (QJ to be adequate. (See also caption of Fig. 

EI). 



In the problem of vortex matter in a trap we encounter 
two specific circumstances, namely a finite-size situation 
and an inhomogeneous density profile. We begin by ex- 
amining the consequences of finite size, by studying the 
system in a cylindrical container with a uniform Pij. 
Such a situation is indeed relevant for the physics of liquid 
^He. Vortex orderings in such a geometry at zero temper- 
ature were studied in JJJ, however the VLL melting for 
this case was addressed only for a planar geometry (see 
e-g. Il3)- Since the melting process in three dimensional 
vortex matter is very different from that in two dimen- 
sions due to the importance of vortex line bending, this 
problem warrants careful consideration. 

Fig. G] shows the results of simulations of vortex mat- 
ter in a cylindrical container with Pij = l^Vij < R. At 
low temperatures, the simulations reproduce orderings 
with circular distortions of VLL near the container wall, 
as predicted for ^He in a zero-temperature treatment of 
the problem 11]. For a large number of vortices the 
system reacquires the hexagonal lattice symmetry away 
from the wall, see Fig. Q] (bottom row). Increasing the 
temperature in the case of small number of vortices (top 
row of Fig. ^ the dominant vortex fluctuations are as- 
sociated with angular displacements, while radially the 
vortex density remains ordered. For a larger number of 
vortices (bottom row of Fig. QJ we find dominance of 
angular fluctuations only for the vortices situated close 
to container wall, while the center of the system does 
not display this phenomenon. The crossover to a uni- 
formly molten vortex system occurs in both cases only 
at a higher temperature. The two-step thermal crossover 
in the vortex pattern we find is analogous to that in two 
dimensions where vortices are point-like objects (see e.g. 
JI3]). There is, however, a principal difference in our 
case, since in three dimensions the VLL melting is ac- 
companied by significant vortex bending fluctuations. 
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FIG. 1: xy positions of vortices in a cylindrical container in- 
tegrated over z-direction, and averaged over every tenth of a 
total of 5 • 10^ MC sweeps. Top and bottom rows have L = 36 
and L = 72, respectively. At T — 0.5, 1.0, we discern circular 
ordering close to the cylinder wall combined with a hexago- 
nally ordered state closer to the center. At T = 1.25, 1.67 we 
observe dominance of angular fluctuations closest to the edge. 

Let us now turn to the case of a harmonic trap (Eq. 



nj. In a uniform and infinite system, the fluctuations can 
cause either VLL melting via a first-order phase transi- 
tion or a second order transition associated with a ther- 
mally induced proliferation of closed vortex loops near 
the critical teraperature where the vortices loose their 
line-tension |3, |l^ • In a finite-size inhomogeneous sys- 
tem, the situation is different. A density gradient in a 
trap may effectively be viewed as a temperature gradient 
in a uniform system. It is clear that for low, but finite, 
temperatures there will be a finite area near the edge of 
the cloud which effectively would be at a high enough 
temperature to feature an annulus of tension-less tangle 
of vortices. This is a phase where the vortex-line ten- 
sion (free energy per unit length) has vanished through 
the proliferation of vortex loops [^ . The boundary where 
this tangle sets in, marks the true boundary of the EEC. 
An issue to be adressed is whether we encounter an ap- 
preciable VL (i.e. tensionful but disordered vortex state) 
layer between the tension- less vortex tangle and the VLL. 
Fig. [2I shows snapshots of vortex configurations in the 
model Eq. Q generated by MC simulations at T = 0.5 
and T = 1.0. For better visualization we choose the 
vortex radius to be 0.4 of the grid spacing so it should 
not be associated with the core size. The sharp bends at 
short length scales result from the presence of a numerical 
grid. Nonetheless, this model has proved to be accurate 
for describing vortex fluctuations at scales larger than 
the grid spacing 



mi: 




FIG. 2: (Color online) Snapshots of vortex configurations in a 
rotating trapped BEG at T = 0.5 (left figures) and at T = 1.0 
(right figures). The top row shows a selection of 16 out of 72 
layers in z direction. The bottom row shows smaller selections 
in the xy plane, but 32 out of 72 layers in z direction. Fluc- 
tuations are minimal in the trap center, and increase towards 
the edge of the trap. A distinct front separating regions of 
ordered and disordered vortices is easily identified. 

We next locate the vortex liquid layer between the 
ordered vortex state and the tension-less vortex tangle 
closest to the edge of the trap. To this end, we take 
snapshots of vortex matter at different temperatures, in 
each case integrating over the z-direction. In a resulting 
image, straight vortex lines will be seen as bright spots 
while bent vortices will be seen as smeared out spots. 
This may be related to experiments, where at least for 



non-equilibrated vortex systems the ^^-integration renders 
vortices essentially indistinguishable |l4'Jl5(. The results 
are given in the upper row of Fig. ISlfor different temper- 
atures. There we can identify regions of rather straight 
and ordered vortex lines and a smeared region. To ob- 



in z direction, defined in a selected region between two 
cylinders of radii R\ and i?2- We do so by applying a 
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FIG. 3: x2/-positions of vortices in a trapped BEG inte- 
grated over z- direct ion. Top row are snapshots, while mid- 
dle and bottom rows are averages of 10^ and 5-10^ MG 
sweeps, respectively. Every tenth configuration has been sam- 
pled. This provides information on the stability of the or- 
dered region and the evolution of the disordered region as 
T varies. Intervortex distance 2ro = 2^/i/MQ and healing 
length ^ = h/y^2Mg2D'n < ca/2, where a is lattice constant 
and c is a constant of order, but less than, unity. Hence, 
2hQ/g2Dn < c^/9, which according to Goddington et al. |l| 
validates our approach. 

tain further insight into the vortex matter in this case, 
we also perform a thermal averaging as in Fig. Q] This 
is shown in the second and third rows in Fig. El By 
averaging over different number of snapshots we identify 
a well-defined boundary between the ordered and disor- 
dered regions. Indeed, in a finite system, averaging will 
eventually produce a complete smearing even in the cen- 
ter of the trap. We observe signatures of this effect in 
the clear differences between the third picture in second 
and third row of Fig. [31 where averaging was made over 
10^ and 5 • 10^ MC sweeps, respectively. Thus, the time 
scale of the fiuctuations in the ordered regions are dra- 
matically larger than those related to the fiuctuations in 
the disordered regions. 

We next investigate the character of the vortex state 
in the disordered region. It is known that in the VL the 
helicity modulus, or equivalently the superfiuid density, 
is zero in any direction [3, |^ • Monitoring of the helicity 
modulus could be employed to identify a region of a pos- 
sible tensionful VL in the above pictures, as explained 
below. In a trapped system the global helicity modulus 
Tz p has no rigorous meaning, due to the non-uniform 
Pij. However, we introduce a modified helicity modulus 



A(r, 



Az if Ri < rij < R2, 
otherwise. 



(2) 



to the model Eq. Q and defining the modified helicity 
modulus as follows, 

t.(i?i,i?2) = ^ l^Pii cos(^,- -Oi- Ai^r^ 

- ^ (E' [^'i si^(^^- -^i- ^'i)]') ■ (3) 

Here, ^ is over all links where A^j is nonzero (depend- 
ing on Ri and R2) and N' is the number of these links. 

In a uniform extended system the proliferation of vor- 
tex loops happens via a second order phase transition 
I3. When loops proliferate, the condensate is destroyed. 
The temperature of vortex-loop proliferation decreases 
with increasing rotation (see Fig. 12 in 8]). Alterna- 
tively, the destruction of the phase-coherence along the 
z-axis is caused by destruction of the lattice order via a 
first order transition. This scenario, if it is realized in a 
trapped system, should be manifest in the shape of the 
helicity modulus as the transition is approached, in that 
it should be significantly different from the case without 
rotation. Namely, one should see a remnant of a first 
order phase transition with a near-discontinuity in the 
helicity modulus, rather than the continuous variation 
characterizing a second order transition driven by a pro- 
liferation of vortex loops. If one were to observe no ap- 
preciable difference in the temperature dependence of the 
helicity with and without rotation, one would conclude 
that the demarcation line seen in the images separates 
an ordered region from tension-less vortex tangle, with 
no discernible VL region. This scenario would imply a 
well defined and regular structure of the boundary of the 
VLL. 

The results for T z{Ri^ R2) are shown in Fig. ^ These 
measurements indeed show that the presence of a rota- 
tion significantly reduces the temperature at which T^ 
vanishes. This reduction relative to the case of no ro- 
tation decreases with increasing i^i,i?2- That is, panel 
(d) is similar to panel (c) (no trap), whereas in panel 
(g) there is little difference between T z{Ri^ R2) with and 
without rotation. Thus, for the latter case the presence 
of vortices essentially does not infiuence T;^, and the de- 
struction of superfiuid density is driven by the prolifera- 
tion of vortex loops. Panel (g) is connected to the three 
leftmost panels in Fig. [2lin the following sense. The dis- 
tinct demarcation line in these leftmost images separates 
a vortex solid from a tension-less vortex tangle with no 
visible tensionful VL region. This is consistent with the 
experiments showing a very regular edffe structure for 
systems with large number of vortices [jj, [3 • Note also 
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FIG. 4: Results for Tz{Ri,R2)- The two top panels show 
thermal depletion of the superfluid density in the model Eq. 
Q (r is the distance from the center of the trap) at the temper- 
atures T = 2.50 (the lowermost curve), T = 2.00, T = 1.67, 
T = 1.25, T = 1.00, T = 0.50. In panels (c)-(g), the 
upper curve (+) is the helicity modulus without rotation, 
while the lower curve ( x ) the helicity modulus with rotation- 
induced vortices with filling fraction / = 1/36 as functions of 
temperature. Panel (c) shows Tz for a cubic uniform sys- 
tem with periodic boundary conditions. The upper curve 
(+) has the properties of a second order transition (helic- 
ity modulus vanishes because of vortex- loop proliferation), 
whereas the lower curve (x) has the finite-size appearance 
of a first order transition (suggesting that the helicity mod- 
ulus vanishes because of vortex lattice melting) 7, 8]. The 
remaining panels (d)-(g) show t40,i?/4), t4i?/4, 2i?/4), 
fz{2R/4:,3R/4), and t^(3i?/4, i?), respectively. The vortex 
plots on the left (obtained at T = 0.50) defines the radii Ri 
and R2 as white circles. Taking parameters from Ref. 2, and 
using Q = (h/ M)Nv / 27: R^ 2] where Nv is the number of vor- 
tices in the trap, we find Q. ^ lOOHz. Since oo± ^ 500Hz [3, 
this puts us well outside the LLL regime. 



the absence of circular distortion for the vortices at the 
edge of the system. On the other hand, by comparing 
panels (c) and (d), we see that T^ and T^ are quite sim- 
ilar. This points toward a possibility of the existence of 
a tensionful VL phase close to the center of the trap. 

We have considered vortex matter in the model Eq. 
Q with uniform and nonuniform density. The uniform 
case features dominant angular vortex fluctuations near 
the wall of the cylinder. Vortex matter in a trapped 
BEG is more complicated, due to density gradients. We 
have identified a number of inhomogenous vortex states. 
Notabley, we find a direct crossover from a vortex solid 
to tension-less vortex tangle with no discernible inter- 
mediate tensionful vortex liquid at low temperatures as 
the edge of the trap is approached. This explains very 
regular edges of VLL at finite temperatures in many ex- 
periments. At higher temperatures, a possible tensionful 
vortex liquid state located between a vortex solid at the 
center and a tension-less vortex tangle closer to the edge 
is identified. Our simulations indicate strong bending 
fluctuations of vortices in this region which may obscure 
its visibility in experiments. An experimental observa- 
tion of the discrepancy between visible vortex numbers 
and rotation frequency (as we predict for disordered vor- 
tex states) has in fact been observed, albeit in an anhar- 
monic trap 15]. 
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